Existence and Vanishing of the Breathing Mode in Strongly Correlated Finite Systems 
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One of the fundamental eigenmodes of finite interacting systems is the mode of uniform radial expansion 
and contraction - the "breathing" mode (BM). Here we show in a general way that this mode exists only under 
special conditions: i) for harmonically trapped systems with interaction potentials of the form 1/r 7 (7 G ffi^o) 
or log(r), or ii) for some systems with special symmetry such as single shell systems forming platonic bodies. 
Deviations from the BM are demonstrated for two examples: clusters interacting with a Lennard- Jones potential 
and parabolically trapped systems with Yukawa repulsion. We also show that vanishing of the BM leads to the 
occurence of multiple monopole oscillations which is of importance for experiments. 
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Strongly correlated finite systems are of high current inter- 
est in many fields. Examples are atomic clusters, e.g., HHIL 
and externally confined systems interacting by repulsive po- 
tentials, such as atomic Bose condensates or strongly interact- 
ing Fermi gases fl3L |4£] , ion crystals in traps |5J, dusty plasma 
crystals, e.g., laLZD, or electrons in quantum dots lElSl]- The 
response of these systems to small external pertubations is 
fully determined by its collective excitations. Among those, 
the radial expansion and contraction - the monopole oscilla- 
tion (MO) - is particularly important since it can be easily 
excited selectively by variation of the confinement |10l] or by 
applying external fields The corresponding frequency 

cjmo can often be precisely measured and may serve as a sensi- 
tive indicator of intrinsic system properties including the form 
of the pair interaction, the trap geometry or the screening 
length and particle charge in complex plasmas 

filill- 
On the theory side, the collective modes are successfully 
analyzed within continuum models, e.g., lfl2i 131. These 
models are applicable to the gas or fluid phases of classical 
or quantum systems where correlations are weak or moder- 
ate. In the strongly correlated crystalline state where the in- 
dividual particle positions rj become separated, such models 
are questionable. Nevertheless, frequently a MO is associ- 
ated with the oscillation of the mean square radius, R 2 (i) = 
Hi r ? (*). OH Q Hi- In three-dimensional (3D) sys- 
tems with Coulomb interaction and harmonic confinement 
with frequency u this MO has the frequency cu M0 — VSujq, 
independent of N, e.g., lfl2l [l4ll . In harmonically confined 
2D systems, also a universal (A-independent) MO was ob- 
served if the interaction is a repulsive power law, ~ 1 jr n (n = 
1,2,.. .), or logarithmic [15]. For other interactions, wmo is 
A-dependent 1U6I1 . 

On the other hand, the response of such a strongly corre- 
lated finite system to small pertubations is fully determined by 
its normal modes. Among those, the mode of uniform radial 
expansion and contraction, the breathing mode (BM), is most 
similar to the MO and both are often used synonymously. 

This common identification of the breathing mode with the 



monopole oscillation assumes that a BM exists always, inde- 
pendently of the confinement and for any pair interaction. In 
this Letter we show, by a general direct analytical investiga- 
tion, that this assumption does not hold for strongly correlated 
finite classical systems. We prove that in fact a BM exists only 
in two classes of systems: i) in parabolically trapped systems 
of any dimension with power law or logarithmic pair interac- 
tion, and ii) in some highly symmetric systems, for arbitrary 
pair interaction. In all other cases, any normal mode devi- 
ates from the BM. As a consequence, then there exist several 
monopole oscillations with different frequencies. 

General existence conditions of the BM. We consider a 
d-dimensional, classical system of A identical particles with 
arbitrary pair interaction v(r) described by the hamiltonian 
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The included confinement potential tj> is isotropic but com- 
pletely general otherwise and may also be equal to zero. We 
assume the existence of a stable configuration (ground state 
or metastable state) r* = (r*, . . . r^) G R dN defined by 
Q2 indicates I 7^ i) 



= ViJ7(r)| P=P 



0V>?l). 



JV 

1=1 



r*u, (2) 



from where the normal modes are excited. For a general nor- 
mal mode analysis we use the eigenvalue equation 



Am? = Tir , 



(3) 



which contains the positive semidefinite Hessian matrix 7i = 
(yiVjU(r)\ T=T *) i j_ 1 N . Accordingly, A = lj 2 > is the 
eigenvalue connected with the mode frequency u>, and f is the 
eigenvector containing the displacement vectors of all parti- 
cles within the mode. 
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To obtain the existence conditions of the BM, we first rear- 
range Eq. (O by evaluating the Hessian matrix H, using the 
isotropy of <fi and the distance dependence of v, which yields, 
for each component i E 1 . . . N 
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(|r«|f'(|r«|)-0'(|r;|)) 



1=1 



(4) 
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Now we insert an eigenvector of the BM given by f oc r* (the 
proportionality constant cancels) and obtain 



A' 



Amr^^'drJIK+^V'drSIK,. (5) 
1=1 

These are the N existence equations of the BM, which just 
have the meaning that, the linearized force on each particle 
from the confinement and the pair interactions has to be purely 
radial and uniform (proportional to r*). To find their solutions 
we separate them into the two conditions. First, radiality is 
equivalent to 



N 



= r* x 



£V(I4I)4, 



(6) 
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which requires Yli v "(\ r ii\) r u = s i r i w i m arbitrary coef- 
ficients Si. These radial equations have two classes of so- 
lutions: Rl) the particle configuration is of special symme- 
try with coefficients Sj determined by the configuration, e.g., 
rotational symmetry with respect to all r*. R2) the pair in- 
teraction v(r) is of specific form. This form can be ob- 
tained by adding a multiple of the tangential components 
of Eqs. (ffjl to Eqs. (6). The result is an equation for v; 
v" (r) = — (7 + l)v'(r)/r, where 7 is an arbitrary real num- 
ber. Its solutions are v(r) ~ l/r 7 [7 7^ 0] and v{r) ~ log(r) 
[corresponding to 7 = 0], and these solutions fulfill Eqs. © 
independently of the particle configuration with coefficients 
given by * = (7 + l)^(|rj|)/|rf|. 

So far we used only the requirement of radiality. To fulfill 
the existence equations <(5j the radial solutions Rl) and R2) 
also have to fulfill the conditions of uniformity which follow 
by multiplying (0 with r* 



Am = 0"(K*|) + 



(7) 



These conditions also have two kinds of solutions: Ul) the 
configuration has another special symmetry, e.g., single shell 
configuration with Si independent of i. U2) in case of the 
radial solution R2) a configuration-independent equation for 



the confinement potential cfr(r) is obtained by using the corre- 
sponding coefficients sf. Am — 4>"{r) + (7 + l)4>'(r)/r. Its 
solutions are 



Am 



2(2 + 7) 
Am 



+ cv(r) 7^-2 (8a) 



Vlog(r) + C7- 2 7 = -2, (8b) 

where 7 is the power of the pair interaction and c is an arbi- 
trary constant. Consequently, there are two basically different 
possibilities for the existence of a BM: 

Universal breathing mode (UBM). The previous analy- 
sis shows that a configuration- and iV-independent univer- 
sal BM exists in the case of harmonically confined systems, 
Eq. ( f8ab , cf)(r) — mw^r 2 /2, with particles interacting via po- 
tentials v(r) proportional to 1 /r 7 or to log(r). For these cases 
the breathing frequency is independent of N and given by 
^bm = V2 + 7 ijJq and wbm = v2 ujq, respectively. For inter- 
action potentials proportional to r 2 , Eq. d8bb . the confinement 
has to be of the form r 2 log(r), where the prefactor determines 
the breathing frequency wbm- In case of particle(s) located in 
the trap origin the confinement is modified containing an addi- 
tional term cv(r), Eq. (8a), or cr 2 , Eq. (8b). These results are 
valid for any real 7 and any dimension and include the result 
of Ref. Il5ll as a special case. 

Further we conclude that no UBM exists for all exponential 
potentials (such as Yukawa, Morse etc.) or non-monotonic 
potentials (e.g., Lennard-Jones). 

Non-universal breathing mode (NUBM). A second class 
of solutions of Eqs. ([5]), the configuration dependent solutions, 
exists in highly symmetric configurations. For instance, 2D 
equally spaced single shell systems or 3D single shell sys- 
tems, which form platonic bodies (with or without particle in 
the center), fulfill Eqs. <(3j for any v and - in these cases 
the system looks the same from every particle's view. Thus, 
for every i £ {1 . . . TV} Eqs. © are identical. Then, since 
no direction is preferred, the sum is proportional to r*, and a 
unique value A > exists fulfilling the existence equations. 
The resulting BM is non-universal, i.e., w BM depends on N, 

Deviations from a BM. We will now verify our predictions 
by a numerical normal mode analysis. To this end, we use the 
ground state configurations of Lennard-Jones and harmoni- 
cally confined Yukawa systems for different N lfl7i [T§] and 
calculate all eigenmodes. For each mode we analyze the de- 
viations from a BM, i.e., deviations from the proportionality 
f = cr*. The first type of deviations arises if fj = ^r-, with 
different Cj for different particles (radial, but non-uniform 
mode). More pronounced deviations exist if there is no pro- 
portionality for all particles, i.e., there exist finite tangential 
velocity components (non-radial mode). Finally, the strongest 
deviations correspond to the case where the mode is not radial 
and contains, in addition, a positive number N ap of anti-phase 
oscillating particles, i.e., these particles move inward while 
the majority moves outward and vice versa. 

We will measure these deviations by computing the number 
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FIG. 1 : QBM of harmonically confined 2D Yukawa systems, a) per- 
fect BM for TV = 40 and nd c — 0. b) purely radial but non-uniform 
QBM for TV = 16 and nd c = 20.0, (<r r = 0, a u = 1.1, jV ap = 0). 
c) non-radial and non-uniform QBM for N = 40 and nd c — 1.99, 
(cr r = 2.77, <7 U = 2.11, 7V ap = 4). The dashed lines mark the shells. 



N- dp whereas radiality and uniformity will be measured by the 
distribution widths normalized to the mean 



max(5i jr — mirier /|#j, r | with Si_ r 



(9a) 



cr u = (ma,xS i:U - niin^^J /5 i)U with S i>u = yy. (9b) 



If, for each mode, cr r and/or a u are non-zero, no BM exists. In 
this case we define the quasi-BM ( QBM) as the mode with the 
smallest of these deviations. Two examples of such a QBM 
are shown in Fig.[T]together with a perfect BM. 

Unconfined 3D Lennard- Jones (LJ) systems. First we 
examine the breathing-type modes of unconfined systems with 
LJ interaction 



v(r) = 4e 



(10) 



which is of relevance, e.g., for molecules and atomic clusters 
id Ufa- While the eigenvalues of the corresponding modes are 
dependent on the chosen energy and length scale, the eigen- 
vectors are not. Thus, a systematic calculation of er r and 
cr u in dependence on N can be performed. The results are 
shown in Fig. [2] where, er r , <r u and iV ap are plotted vs. N 
in the range from 3 to 150 particles. One clearly sees that 
most of the systems show substantial deviations from radi- 
ality and uniformity. The smallest deviations are observed 
for N = 55, 135, 147, which are highly symmetric multi- 
ple shell systems with full icosahedral symmetry lfl7tl (see 
dashed-dotted lines): the values are er r = 0, iV a p = and 
<7 U = 0.2, 0.24, 0.23, respectively. Thus we confirm that LJ 
systems do not possess a UBM, and a NUBM exists only for 
N = 3,4, 6, 13 (see dashed lines in Hig. [2j and the trivial case 
N = 2. These particle numbers correspond to single-shell 
configurations, which are either planar (N = 3) or form pla- 
tonic bodies: tetrahedron (N = 4), octahedron (N — 6) and 
icosahedron with a central particle (N = 13). Other platonic 
configurations do not show up. 

Parabolically confined Yukawa systems. As a second ex- 
ample of broad practical interest, e.g., for spherical dusty plas- 
mas and colloidal systems, we consider the breathing-type 
modes in a parabolically confined one-component Yukawa 
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FIG. 2: QBM of unconfined LJ vs. particle number, a) Deviations 
<j r from radiality. b) Deviations a u from uniformity, c) Number N ap 
of anti-phase oscillating particles. Cases of NUBM are marked by 
dashed lines (see also Fig. [4]l and the QBM of the icosahedral sym- 
metric systems by dashed-dotted lines. 
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using as length scale d c — (2Q 2 /mw^) 1 ' 3 - the stable dis- 
tance between two charged particles in the absence of screen- 
ing tM- F° r K ~ * (Coulomb system), a universal breathing 
mode is observed, in agreement with the previous analysis, 
see Fig.[T]a. In contrast, for n > our simulations confirm 
that no UBM exists [19]. Two typical examples of the QBM 
of 2D Yukawa systems are presented in Fig.[T]clearly showing 
deviations from radial and uniform motion. While in Fig.[T]b 
still a purely radial motion is observed, in Fig.Q]c, the mode 
closest to the BM shows striking deviations, including several 
particles moving in tangential direction. 

Let us now analyze the behavior of the QBM as a func- 
tion of the screening length. Surprisingly, small variation of k 
leads to an irregular behavior of <r r and/or a u with very sharp 
peaks, cf. Fig. [3] showing the properties of the QBM for the 
example N = 40 in 2D and < nd c < 2. The origin of 
the peaks is easily understood: they arise whenever the QBM 
switches from one eigenmode to another. I.e., one mode hav- 
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FIG. 3: (Color online) Deviations a, (solid line) and a u (dashed line) 
of the QBM for harmonically confined 2D Yukawa systems with 
N — 40 vs. screening strength Kd c . The inset shows <j r and a a 
of modes 46 (red line, numbers correspond to the mode frequencies 
in descending order) and 47 (blue line), which have the smallest de- 
viations. The QBM is presented by dots. 
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FIG. 4: Frequency of all existing NUBM of 3D LJ (right axis) and 
2D and 3D Yukawa (YU) systems with Kd c = 1 (left axis). 



ing the smallest deviations from the BM in a certain k range 
loses this role to another mode at some critical value of the 
screening parameter. This can be seen in the inset of Fig. [3] 
showing the QBM peak in the region 1.965 < Kdc < 2. 

Furthermore, we consider the NUBM of LJ and Yukawa 
systems more in detail. All existing BM (for a single value 
of k) due to highly symmetric configurations are shown in 
Fig- S] Note that, in the case of harmonically trapped Yukawa 
systems, the values of N exhibiting a NUBM depend on n 
and on the dimensionality. Further, the simulations confirm 
that the frequencies are non-universal, i.e., depend on N. 

Finally, let us now return to the monopole oscillation dis- 
cussed in the introduction. Considering small harmonic par- 
ticle oscillations with eigenfrequency u> around the stationary 
positions the mean square radius becomes R 2 (t) = Rq + 
2sinwi r* f + (9(f 2 ). If the BM exists, the d-N— dimensional 
vectors f bm and r* are parallel and there exists a MO with 
ujmo = <^bm, which is unique, owing to the orthogonality of 
all eigenvectors. In contrast, if no BM exists, generally many 
eigenvectors f ^ with frequency ui^ will have a non-vanishing 
projection on r*, each giving rise to a MO with different fre- 
quencies cli^o — w fc an d amplitudes. In this case, an exter- 
nal perturbation of the system which excites radial particle 



motions will excite all those monopole oscillations at once. 
The resulting complex spectrum is a fingerprint of the internal 
properties of the system, in particular the form of the pair in- 
teraction and the confinement, and is expected to be a sensitive 
experimental diagnostic of strongly correlated systems. 

In summary, we have derived the general existence con- 
ditions for the BM in strongly correlated finite classical sys- 
tems of any dimension with arbitrary confinement and inter- 
action potential. The existence of a BM is restricted to a small 
class of systems: i) parabolically confined systems with a 
power law or logarithmic interaction potential, or ii) partic- 
ular, highly symmetric configurations, independently of the 
interaction. Systems with exponential or non-monotonic po- 
tentials (including Yukawa, Morse or Lennard-Jones systems) 
not having such a special configuration do not possess a BM. 
In these systems all normal modes exhibit non-uniform and/or 
non-radial motion which has been illustrated for several rep- 
resentative examples. 
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